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Abstract 



We propose a novel method to measure the fractal dimension of a submono- 
layer metal adatom system grown under conditions of limited diffusivity on 
a surface. The method is based on measuring the specular peak attenuation 
of He atoms scattered from the surface, as a function of incidence energy. 
The (Minkowski) fractal dimension thus obtained is that of contours of con- 
stant electron density of the adatom system. Simulation results are presented, 
based on experimental data. A coverage dependent fractal dimension is found 
from a two-decade wide scaling regime. 
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I. INTRODUCTION 



Fractal dimension (FD) [l|] has proved to be a very useful concept for the description of 
certain disordered systems 0,^. The experimental determination of this quantity is therefore 
an important challenge. Scanning tunneling microscope (STM) studies of adsorption and 
homoepitaxial growth of metal atoms on metal surfaces have recently revealed the formation 
of fractal-like structures, under conditions of limited adatom diffusivity and low deposition 
flux The characterization of the structures formed during these growth processes is 

a major objective of surface science, and the determination of FD plays an important part 
in meeting this objective. At present, STM is virtually the only available method for mea- 
suring the FD of such 2D systems. This involves image analysis of the STM topograph, as 
performed, e.g., by Hwang et al. for Au/Ru(0001) and Brune et al. [H for Ag/Pt(lll). 
The STM method has, however, several experimental limitations. In particular, the method 
is inherently local, and the statistics for characterizing disordered surfaces, including fractal 
structures, may not be the best. In general, to probe these systems at the microscopic level 
requires resolution of several A. Scattering techniques used conventionally for the measure- 
ment of fractal dimension at such high resolutions, such as X-rays and neutron diffraction 
1^, are not applicable, due to their penetration into the bulk. An attractive alternative is 
presented by He scattering, which is particularly well suited to the study of disordered sur- 
faces [H as it is not destructive, probes solely the surface, and provides excellent statistics: 
A single beam can cover an entire surface of close to a cm^ at once, compared to merely 
several (lO^A)^ with STM. 

In the present Letter, we propose a novel method (SIMF - Specular Intensity Measure- 
ment of Fractality) by which a He scattering experiment can be used to measure the FD 
of submonolayer disordered surfaces. In sharp contrast to STM, which due to its inherent 
local nature observes individual clusters, the proposed method measures the FD of an en- 
tire disordered surface. Thus the fractal property observed by a SIMF experiment, differs 
significantly from the STM one. Another difference concerns the physical interpretation of 
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what is actually measured: STM probes the electron density near the Fermi level. In the 
repulsive regime, He experiences an exchange interaction with other electrons: It is scat- 
tered from regions that represent typically the low density tail of the surface electrons. In a 
semi-classical picture, at a given incidence energy. He is scattered from contours of constant 
potential. Using the Esbjerg-Norskov-Lang relation 

Vrepix,y,z) = ap{x,y,z) (1) 

for the repulsive part of a He/metal interaction, where a is a system-dependent constant and 
p{x, y, z) the surface electron density, it follows that it is the FD of constant electron density 
contours which is measured in a SIMF experiment. These two differences illustrate that 
SIMF, which offers a global view at the electron density tail level, is in fact complementary 
to the measurement of FD by STM, which measures locally the electron density at the Fermi 
level. More explicitly, since SIMF looks at contours of constant and low electron density, it 
observes an object which is quite different from the underlying geometrical set of adatom 
positions. This object is the "thickened set", or Minkowski Sausage" of the adatoms 
(see inset of Fig.^. SIMF is simple to implement experimentally: All it requires is the 
measurement of specular intensities, over a wide range of incidence energies. This has to 
be performed for calibration on a very low coverage system (ideally a single adatom), and 
then on the system of interest itself. A straightforward analysis subsequently yields the FD. 
The method is based on the energy-dependence of the scattering cross section (CS) of He 
scattering by an adatom P,|10|- This dependence can be exploited to observe the adatom 



system at different resolutions, thus providing the variation of yardstick sizes necessary for 
FD measurements. 

II. THEORY 

Consider first the measurement of the CS of a single adatom adsorbed on an otherwise 
flat surface (as is the case to a good approximation for Pt(lll)). Assume also that the 



3 



surface is perfectly rigid. Any off-specular intensity must be attributed to scattering by the 
adatom ||^ , which allows for a convenient measurement of the CS as 

/(/t,) = /oe-"'^=('=^). (2) 

Here / and Iq are, respectively, the specular intensities for scattering from a surface with a 
single adatom and a flat surface; n is the number of sites per unit-cell; 6 is the coverage; 
is the adatom CS; and kz is the incoming He wavenumber, assumed for simplicity to strike 
the surface at normal incidence. One can accordingly associate a scattering radius, Rs{kz) = 
{T,s{kz) /TrY^"^ with a single adatom. (Formally Bragg scattering must be counted along with 
specular scattering, but this effect is negligible for the fiat metal surfaces considered here. 



We have shown elsewhere [jTO[ that Eq.(g) is rigorously equivalent to the formal definition 
of the scattering cross section of a single adatom). Next consider the CS measurement of a 
disordered submonolayer of adatoms. Suppose that every collision with an adatom results 
in off-specular scattering. This may be the case if the coverage is sufficiently low, so that no 
large compact regions have formed; or if the ad-islands are "dendritic" (fractal), so that the 
width of their "arms" is only several atoms across; or if the adatoms are non- metallic (e.g., 
CO), so that they present a highly corrugated surface to the incident He beam. Then the 
probability of hitting the adatom system is just the off-specular intensity. But this is also 
clearly the fraction of surface occupied by the adatom system. Thus 

Asikz) = l-Iikz), (3) 

where / is now the specular intensity for scattering from the adatom system, and Ag is the 
total CS of the diffusely scattering adatoms, normalized to the total available surface area. 
Thus As cannot be considered the total CS of large compact metallic islands, which scatter 
diffusely only from their rim. This should, however, not be regarded as a limitation of SIMF: 
The method is designed to be applied foremost to fractal structures, and these are k;nown 
to have "arms" that rarely exceed 4-5 atoms across If specular scattering from the 

interior of ad-islands is important, the adatom system is likely to contain large compact 



regions and should probably be viewed as a "fat fractal" |TT|, i.e., a fractal with non-zero 

i.e., 



Lesbegue measure. Such fractals are characterized by their "exterior dimension" |T2| 
the fractal dimension of the rim of the islands. But since it is the rim which accounts for 
the off-specular intensity, one can still apply SIMF, and find that it measures to a good 
approximation just this exterior dimension. Thus the method is seen to naturally provide 
the proper characterization of fat fractals as well. A point of concern regarding Eq.(|^) is that 
it implicitly employs the "geometric overlap principle" of adatom CSs. This has been 
criticized for adatom pairs |]13|. For other systems |l8| Jl^ , [T5[| it is strongly corroborated by test 



cases. We accepted it here, by assumption, as a method valid for all systems, for sufficiently 
low coverages. For some systems, at high coverage, systematic higher-order corrections to 
Eq.(^ may be necessary, which in our belief can be introduced. However, in this study we 
shall not deal with these complications. 

III. RESULTS AND DISCUSSION 

We proceed to present the results of our analysis of combined experimental and simulation 
data. To simulate (presently unavailable) experimental data, we performed He scattering 
calculations within the Sudden approximation (SA) |16| on configurations of Ag/Pt(lll) 



generated by kinetic Monte Carlo (KMC) simulations IT^, in accordance with the Deposi- 
tion, Diffusion and Aggregation (DDA) scheme [|18|. This scheme has been very successful in 
modeling thin-film-deposition experiments. A typical configuration generated by our simula- 
tions at a relatively high coverage of 27% is shown in the inset of Fig.|, and is reminiscent of 
the STM topographs of Refs. PJ^. Experimental CS data was available for scattering from a 



single Ag/Pt(lll) adatom [jT9[. The He/Ag/Pt potential (to be published separately in Ref. 
poll), unavailable at present from other sources, was obtained by fitting SA-CS calculations 
to this data. (To test its reliability, the same potential was used independently to reproduce 
an experimental angular intensity distribution for scattering from a system of polydispersed, 
small Ag/Pt(lll) clusters, which it did with impressive accuracy ||2^). The fit is shown in 
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Fig.|T]. Clearly, this potential is approximate only. However, we believe it is realistic enough 
for the purpose of simulating scattering data. In any case, the method does not depend on 
any specific potential: The ultimate test of SIMF will have to be experimental. The single- 
adatom CS, ^^(/c^), shown in Fig.|I], was extracted from the specular intensity according to 
Eq.(^, yielding Rs{kz). As seen, beyond fc!" ^ lObohr"^ lO'^K, T,s{kz) tends to a minimal 
limiting value S^ = 7r/2^, with /2m = 4.55A. The Pt(lll) lattice-constant is 2.77 thus the 
CSs of nearest-neighbor Ag adatoms always overlap, and the set appears as a continuum to 
the incoming He. This is shown in the inset of Fig.|, where a circle of radius Rm is drawn 
around each nucleus, and the resulting contour is shown. This is an approximate view of the 
contour of constant electron density observed by He atom incident at > A;™, and of the 
fractal whose dimension is measured by SIMF. The adatom system CS, As{kz), was found 
according to Eq.(|) from the specular intensities calculated at various coverages, shown in 
Fig.|. 

Once the functions Rs{kz) and As{kz) are known, the FD can be found. The algorithm 
we propose is based on the Minkowski |^ definition. The Minkowski cover M{R) of a set X 
is obtained by placing circles of radius R centered at every point in X. One then calculates 
the area A{R) of the union of all circles. The FD is found in practice as the slope at small 
R values of the log-log plot 



which presumes that a constant slope (i.e., scaling) is observed over at least one decade ^1 . 
In the present context, a natural choice of the set X is the union of adatom CSs at the 
high /cz/minimal radius (-Rm) limit, shown in the inset of Fig.^. Since in this limit the CS 
is determined by the repulsive part of the potential, semi-classically, the CS circumference 
corresponds, through Eq.(0), to a contour of constant electron density. Therefore, this 
choice yields the FD corresponding to such contours of the Ag system, where a high-energy 
He atom would be stopped. This is an inherent property of the adatom system, uniquely 
accessible to He. Other probes (e.g., STM) would measure the FD of another feature of the 



adatom system. The remarkable property of He scattering, is that it is capable of providing 
the Minkowski cover of the adatom set, through a variation of the incidence energy, which 
is equivalent to a change of the scattering-radius Rg. Thus, the process summarized by 
Eq.(^ can be carried out in an He scattering experiment in two steps: First, the CS of an 
individual Ag adatom is measured over a wide range of incidence wavenumbers, at least as 
far as the saturation point A;™. This then yields Rs{kz), with the minimal radius measured as 
Rm = Rsik"^)- Second, Asikz) is measured under the same conditions of variable incidence 
wavenumbers. Eliminating between Ag and Rs = Rm+R yields the scattering-determined 
area of the Minkowski-cover, As{Rs)- From here the calculation of the FD proceeds according 
to Eq.(|). (Note that it is the "Minkowski-radius" , R, and not the full cross- sectional radius 
Rg, which appears in Eq.(^)). We applied this procedure, for step 1, to the CS fitted to the 
experimental Ag/Pt(lll) data (Fig.|l|), and for step 2, to our scattering simulation data at 
various coverages. The results are shown in Fig.^ It is noteworthy that 2 decades of linear, 
scaling behavior are observed. The existence of a large scaling regime is an "empirical" 
proof of the fractality of our systems, and the ability of He scattering to measure it. It 
is tempting to compare the scattering results to a simple geometrical model, where the 
Minkowski procedure is applied to circles of radius Rm, centered at the adatom positions. 
However, such a naive approach is misleading: The quantum CS has no simple geometrical 
counterpart. The CS is "fuzzy", i.e., a He atom is scattered off-specularly from within the 
CS only with a probability. Work incorporating a distribution of geometrical radii, allowing 
a comparison with the CS, is currently in progress. 

The corresponding slopes, or FDs, are displayed as a function of coverage in Fig.|[ The 
FD rises monotonically with coverage and is close to the embedding-space dimension of 
2. This high FD is expected, if the large value of the individual CS (Fig.|l|) is taken into 
account: As seen in the inset of Fig.^ due to this large CS the contour has a rather smooth 
structure (compared, e.g., to a typical, highly ramified, DLA cluster), thus enclosing an 
almost "ordinary" 2D set. The non-universality of the FD (i.e., its dependence on coverage) 
is a definite prediction of the theory, and differs significantly from the well-known result for 
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individual DLA clusters, the only fractal systems observable so far by STM , whose FD? 



1.7 ||2^ is universal. However, it does agree with the results obtained by Jensen et al. |T8| for 
the DDA model: There too, a strong dependence of the FD on coverage is observed. More 
importantly, as previously mentioned, Rm > cl, the Pt(lll) lattice constant. This implies 
that the geometry probed by the He beam is much less porous that that of DLA clusters. 
Gaps in the latter can easily be smoothed out by overlapping CSs, and disjoint DLA clusters 
may in fact become one from the He point of view. The rise in FD can be understood as an 
increasingly larger portion of the surface being covered by overlapping CSs, as the coverage is 
increased. Taken together, these observations should explain why the geometry of constant 
electron density contours, probed by a He beam, has a very different fractal nature from 
that of DLA clusters. We expect actual He scattering experiments to produce similar result 
to Figs.|I|-|. 



IV. CONCLUDING REMARKS 

In summary, a simple new procedure for the measurement of FD of disordered sub- 
monolayers by He scattering has been proposed. First the specular intensity as a function of 
incidence wavenumber for a single adsorbate must be measured, from which the scattering- 
radius Rg is determined according to Eq.(|]). Then the specular intensity is measured over 
the same incidence wavenumber range for the adatom system of interest, and its CS, Ag, is 
found according to Eq.(|^). Eliminating the wavenumber between these two quantities, the 
function As{Rs) is plotted on a log-log scale according to Eq.(^). The slope of this plot is 
the FD. If the adatom system contains large compact regions, the calculated FD should be 
interpreted as the exterior dimension of a fat fractal. 

SIMF measures the FD of constant electron density contours of an adsorbed metal sys- 
tem. This observation raises some interesting questions for future research, which is beyond 
the He scattering context: For example, what is the effect of this fractality on bulk proper- 
ties, such as surface conductivity and elasticity? 
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Finally, the present article does not deal with the effect of adatoms adsorbed near step 
edges. These adatoms have different CSs from those adsorbed on terraces, and hence require 
a separate calibration of ^^(/c^). This issue is currently under investigation. An initial 
experimental effort to implement SIMF for Ni/Cu(lll) is under way [ p3| . 
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FIGURES 

FIG. 1. Cross section of single Ag atom on Pt(lll) as a function of incidence wavenumber. 



Circles are experimental results |19|, accurate to within 20%. Dots connected by solid line are 
calculated CS values, using an LJ 6-12 potential with parameters fitted to the experimental data. 
The CS saturates at k"^ ~ 10 bohr^^. Small interference oscillations due to particles striking the 
flat surface and the adatom tops can be observed. These measurements and simulations constitute 
step 1 in the SIMF algorithm. 

FIG. 2. Specular intensity of He scattered by fractal adatom islands of Ag/Pt(lll), for cover- 
ages of 3%-27%, averaged over 30 configurations at each coverage. Interference oscillations increase 
in magnitude with the coverage. The specular intensity overall increases with wavenumber, reflect- 
ing the decrease in total cross section of the Ag islands. These simulations constitute step 2 in the 
SIMF algorithm. 

FIG. 3. Plots of log[As{Rs)]/ R'^ vs log{Rm/ R), for various coverages, used to determine the 
FD (note that the oscillations from Figs.^,^ have cancelled out). Two decades of scaling behavior 
are observed, indicating that He is indeed capable of measuring the FD of the scattering object. 

FIG. 4. He scattering prediction of fractal dimension as a function of coverage, obtained from 
the log-log plots in Fig. 2. Error bars are due to linear regression. Inset: Dots are nuclei positions 
in typical configuration obtained by KMC simulations on a Pt(lll) hexagonal surface (periodic 
boundary conditions apply). Coverage is 27%. Contour represents view of electron density by He 
atom incident with > k^. That is, each adatom is drawn to scale with a cross-sectional radius 
of Rm = 4.55A. According to SIMF, this fractal has a dimension of 1.98. 
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